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rj , Classical particle mechanics on curved spaces is related to the flow of ideal fluids, by a dual 
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Abstract 



jj-^ ' interpretation of the Hamilton-Jacobi equation. As in second quantization, the procedure 

^~>' relates the description of a system with a finite number of degrees of freedom to one with 

inflnitely many degrees of freedom. In some two-dimensional fluid mechanics models a duality 
transformation between the velocity potential and the stream function can be performed 
relating sources and sinks in one model to vortices in the other. The particle mechanics 
counterpart of the dual theory is reconstructed. In the quantum theory the strength of sources 
. , and sinks, as well as vorticity are quantized; for the duality between theories to be preserved 

^— V ' these quantization conditions must be related. 
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1 Particles 

The free motion of a classical particle with unit mass, moving in a smooth space 
with metric gij{x) is described by the Lagrangean 

L = -gij{x)x''x\ (1) 

where as usual the overdot represents a time-derivative. The Euler-Lagrange 
equations imply that the particle moves on a geodesic: 

^=x' + T„:±^i^ = 0. (2) 

The canonical formulation of this theory is constructed in terms of the momenta 

and the hamiltonian 

H = ^g'^p,p,. (4) 

The time-development of any scalar function F{x,p) of the phase-space co- 
ordinates is then determined by the Poisson brackets 

dt ' 5x* dpi dpi dx'- ' 

In particular the Hamilton equations themselves read 

x' = — ■■ = -^ (6) 

dpi ' * 5x* 

A third formulation of the classical theory is provided by Hamilton's principal 
functionFI S{x,t), which is the solution of the partial differential equation 

— =-H{x,p = VS). (7) 

For the case at hand this Hamilton-Jacobi equation takes the form 

^=--9'^V.SV,S. (8) 

Particular solutions S are provided by the action for classical paths x*(t) obey- 
ing the Euler-Lagrange equation @), starting at time r = at an initial point 
x*(0), and reaching the point x*(t) = x* at time t = t: 



S{x,t) = / dTL{x,x) 
Jo 



(9) 



^The terminology follows ref. [ 



An example of the class of theories of this type is that of a particle moving on 
the surface of the unit sphere, S^. A convenient co-ordinate system is provided 
by the polar angles {6, ip), in terms of which 



L 



(9,^) = ^ [e^ + sm^9ip^ 



for a particle of unit mass. The corresponding hamiltonian is 



2 , Pf 

sm ( 



2 ' 



with the momenta and velocities related by 



Pe 



p^ = sm a ip. 



(10) 



(11) 



(12) 



The second equality (|ll]) relates the hamiltonian to the Casimir invariant of 
angular momentum, the components of which are constants of motion given by 

Jx = -sin (fpg- cos ipctg 9 p:^, Jy = cos ippg - sin fctg 9 p^, Jz = P^- 

(13) 
The geodesies on the sphere are the great circles; they can be parametrized by 

cos0(t) = sinasinw(r — r^,), tg ((/^(r) — 99,,) = cosa tgu;(r — r^,), (14) 

where a is a constant, and r* and ip^, are the time and longitude at which the 
orbit crosses the equator: 9^ = 7r/2. On these orbits the angular frequency is 
related to the total angular momentum by 



UJ 



2H = y, 



(15) 



Observe that, for an orbit reaching the point with co-ordinates (0,<y9) at time 
T^, +t, the following relations hold: 



cos a; = sin9 cos{ip — (p^), sinuit = Jl — sin 9 cos2(c/5 — (p^), 



cos 9 



(16) 



sma 



1 — sin^ 9 cos^(c/3 — (f*) 



The last equation implicitly describes the orbit 9{ip), defining a great circle 
which cuts the equator at 9 = 9^ = 7r/2 and ip = ip^, at an angle a defined by 
the direction of the angular momentum: 



J. 



cos a. 



J± 



sin a, J± 



VJ2 VJ2 

The Hamilton-Jacobi equation for this system reads 



'^x ' '^y 



dS_ 
'dt 



dS\' 



+ 



1 



dS 



89 J ' sin^9 \dip) 



(17) 



(18) 



The solution corresponding to the orbit (14) is 

1 



S{6, (p,t) = — arccos [sin 6 cos{(p — 99*)] , (19) 

^0 



ich satisfies the equations 


dS 


(jcosO cos(ip — (/?*) 


de 


/ 9 

yl - siir6cos^{ip - ip^) 


as 


wsin^ sin(c^ — ip^) 


dip 


"^ / 9 ' 

y 1 - sin^ 9 cos'^{ip - ip^) 


dS 
dt 


2 



(20) 



In this approach, the expressions on the right-hand side are obtained by defining 
uj via the last expression, in agreement with (p!6|). The same principle of energy 
conservation/time-translation invariance implies that S does not depend on t^,. 

2 Fluids 

The Hamilton-Jacobi equation (^) can itself be obtained in a straightforward 
way from a variational principle: introduce a Lagrange multiplier field p{x) and 
construct the action functional 

A{p,S) = J dt J d'^x^p^^ + ^g^^ViSV,Sy (21) 

The square root of the (time-independent) background metric has been included 
to make the integration measure invariant under reparametrizations. Of course, 
we could absorb it in the definition of Lagrange multiplier field, but then p would 
transform as a density rather than as scalar. 

The Hamilton-Jacobi equation follows by requiring the action to be station- 
ary w.r.t. variations of p: 



1_6A _dS_ 1 
^J^~~dt^2 



— — = — + o 9*' V,5 V,-5 = 0. (22) 



On the other hand, the stationarity of A{p, S) w.r.t. S gives 

+ V, {g^'pVjS) = 0. (23) 



J_SA _dp 

^Js~~di 



This equation can be interpreted as the covariant equation of continuity for a 
fluidg with density p and local velocity 

Vi = ViS => ^ + Vi [pv') = 0. (24) 



^For background, see e.g. ref.| 



In this interpretation the gradient of the Hamilton-Jacobi equation gives the 
covariant Euler equation 

^ + v^V,v. = 0, ^jv, = ^-r^^v,. (25) 



Eq.(p4|) states that the fluid flow is of the potential type. Indeed, in the absence 
of torsion the Rien 
vorticity vanishes: 



of torsion the Riemann-Christoffel connection F- • '^ is symmetric and the local 



V^Vj - Vjv, = 0. (26) 

For the fluid flow to be incompressible, the velocity field must be divergence 
free: 

V ■v = AS = 0, (27) 

where A = g^^ '^■i^ j is the covariant laplacean on scalar functions over the 
space. It follows that the number of incompressible modes of flow on the mani- 
fold equals the number of zero-modes of the scalar laplacean. For example, for 
flow on the sphere S"^ (or any other compact Riemann surface) there is only 
one incompressible mode, the trivial one with u* = everywhere. 

For a given geometry gij{x), the solution of the Hamilton-Jacobi equation 



(^), ( [2^ ) provides a special solution of the Euler equation (25); for a conservative 
system: dS/dt = —H = constant, it implies dvi/dt = and v^VjVi = 0. Ac- 
cordingly, this solution describes geodesic flow starting from the point {9^, c/?^,). 
To turn this into a complete solution of the fluid-dynamical equations (^) , 
(P5|) it remains to solve for the density p. The equation of continuity takes the 
form 

^+V,{pV'S) = 0. (28) 

It follows that a stationary flow, with p not explicitly depending on time t, is 
possible if 

V • (pVS) = 0. (29) 

In addition to the trivial solution p = po = constant, v = VS/m = 0, it 
is possible to find non-trivial solutions of equation ( p9[ ) for spatially varying 
density p. As an example, we consider flow in a 2-dimensional space; in this 
case one can introduce a generalized stream function T{x,t), dual to the fluid 
momentum, and write 



1 

Then for theories of the type (pi) : 



pV'S = -^ e'NjT. (30) 



" ^{VSr 2H^g ■ ^^ ^ 

With H constant, the factor 2H in the denominator can be absorbed into the 
definition of T = T/2H, and hence the density is given by 

P=^ e'^ViSVjf = -^ e'^ViVjf, (32) 



for the pseudo-scalar function T the gradient of which is dual to pSJ S. Note 
also, that eq.® implies VS" • VT = t; • VT = 0. 

As an illustration, we again consider the unit sphere S'^ . The velocity field is 
given by the momenta (^) per unit mass: 

uJcosO cosiip — iftf) usmO smiip — ip^,) 

V0 = , =, % = . ^. (33) 



1 — sin 9 cos"^ {ip — if ^) yl — sin 6cos'^{(p — ip^) 

Taking into account that on the sphere the non-vanishing components of the 
connection are 

V = i^' rj = -sin0cos0, (34) 

a straightforward calculation shows that indeed 

(ft) ' 

VjV^ = CJ^ v^VjVi = 0, TT " °' ^^^^ 

The first two equations actually imply v^^ViVj — 'VjVi) = 0, in agreement with 
the absence of local circulation. From these results it follows, that the flowlines 
are geodesies (great circles) given by eq. ([l6|) , and stationary. 

For the gradient of the stream function T to be orthogonal to the velocity 
field (|33|), it must satisfy the linear differential equation 

vVT = <^ tg{ip - ip^)V^T = sine cose VeT. (36) 

The general solution can be obtained by separation of variables, and is a func- 
tion of the single variable: T{e, ip) = f (y), with y = tge sin(c/? — (/?*)= ctga. 
For such a scalar field 



s'm{ip — (/?*) 
cos^ 



v.r = ^£^ /'(y)|,=.,. > v,r = tge cos(^ - ^.) /'(y)|,_,^„ 



(37) 



The corresponding density p is then 



p{0,ip) = M = ^i /'(2/)|,= .,.- (38) 

cos e to sm a cos e >y ctgu 

The simplest, most regular solution is obtained for p{a) = p* sin a: 

p^ sin a p^ 

P(G, V) = T- = I (39) 

^°^^ ^Jl-smHcos\^-^.) 

This solution corresponds to 

T{e,p) =ujp^a{e,p) 4^ f{y) =ujp^arcctgy. (40) 

Observe, that in this case T, like a, is an angular variable; indeed, a increases 
by 27rn on any loop winding around the point (0 = it/2; ip = (p^,) n times. 



The solution ( |39D possesses singular points at = tt/2, ip = Lp^ + mr, 
corresponding to a source for n = 0, and a sink for n = 1. This can be 
established from the expression for V • w: 

ujsm9cos{(p-(p^) 
v ■ V = , (41) 



1 — sin 6 cos^ (c/9 — 9?^, ) 

which becomes (+cxd, — cxd) at the singular points. However, a more elegant way 
to establish the result, is to make use of the stream function ( ^0| ) and consider 
the flux integral 

^T) = jdlpvn, (42) 

representing the total flow of material across the closed curve V per unit of time. 
Consider a contour V winding once around the singularity at {6 = 7r/2; ip = ip,,)] 
on such a curve a increases from to 27r. Then 

$(r) = / ^Eijpv'dx^ = i ViTdx' = 2-KUjp^. (43) 

This represents the total flow of matter from the hemisphere centered on the 
source at {9 = 7r/2; p = ip^) to the hemisphere centered on its antipodal point, 
the sink at {6 = vr/2; 99 = (/?*+ vr). 

3 Vortices 

The dual relationship between the velocity potential S and the stream function 
T suggests to study the dynamics of a fluid for which T is the velocity potential: 

v^ = - ViT. (44) 

P* 

The constant /?,, has been included for dimensional reasons. Like before, this 
velocity field is stationary: dvi/dt = 0, but it is not geodesic. Indeed, the 
velocity field describes motion under the influence of an external potential; 
specifically: 

v.Vvi = lViV^ = ^ ViiVTf = -^2 V^{pVSf. (45) 

2 2pi 2pi 

Here p{x) and S{x) denote the previously defined functions mapping the man- 
ifold to the real numbers — e.g. (19) and (p^) for fluid motion on a sphere — 



irrespective of their physical interpretation. Now again, as (VS") = 2H = uj 
constant, it follows that 



u^ 



2pi 
Combining eqs.(p5|) and 



^■^'^i = Wj2^iP = -^i^- (46) 



2„2 



l-' = -ih-ho) = ^, (47) 



where h represents the external potential. Because of the potential nature of 
the flow, eq. (^^, the local vorticity again vanishes: ViVj — VjVi = 0, but as 
eq.(|43|) shows, this is not necessarily true globally. Indeed, in singular points 
of the original geodesic fluid flow (with sources/sinks), the dual flow generally 
has vortices/anti- vortices. 

Continuing our example from the previous sections, we can illustrate these 
results in terms of flow on the unit sphere, for which T/p^ = uja and Vi = ujVia: 



uism{ip — (f^) ujsm6cos6cos{ip — ip:t:) 

1 — sin^ 9 cos'^{(p — if^) ' 1 — sin^ 9 cos^{ip — 99^,) 



It follows, as expected, that 

v' = -\yo.f = i^ = . "' -. (49) 

p% 1 — sm 9 cos^[(p — (ft,) 

A further remarkable property, is that the dual flow is divergence free: 

V • u = -^ Aa = 0, (50) 

where-ever v is well-defined; obviously, the result can only be true because of 
the two singular points {9 = 7r/2;93 = (p^,) and {9 = 7r/2;c/3 = 93* + vr), where 
Vi and its divergence are not well-defined, i.e. topologically the velocity field is 
defined on a cylinder, rather than a sphere. These two points are centers of 
vorticity, as follows directly from eq. (|43|) , which in the present context can be 
rewritten as 

® Vidx' = 2-Kuj, (51) 

for any closed curve T winding once around the singular point ('7r/2, ip^); as this 
curve also winds once around the other singular point in the opposite direction, 
they clearly define a pair of vortices of equal but opposite magnitude. 

As the flow is divergence free, it follows that in this case there can be non- 
trivial incompressible and stationary flow modes: for constant density pi one 
has 

^ = 0, Vpi = => V-(pi?;) = 0, (52) 

and the equation of continuity is satisfied. 

The nature of the flow lines deflned by eq.(p8|) is clear: they are parallel 
circles of equidistant points around the centers of vorticity. On these circles 
the velocity is constant in magnitude, implying by ( |49| ) that sin cos (99 — 99*) = 
cos/? = constant. For example, for 99* = we get x = cos/? = constant; the 
flow line then is the circle where this plane of constant x cuts the unit sphere. 
On these flow lines 

vq = —ujism[p — ipt,), Vip = —iOicos l3cos9, (53) 



with 



v"^ CO u 



to 1 — sin^ 9cos'^{(p — (p^:) sin^ /? 



4 The dual particle model 



Having clarified the nature of the (incompressible) flow described by the dual ve- 
locity potential T/ p^, , we now reconstruct the corresponding particle-mechanics 
model for which T/p^: is Hamilton's principal function. From eqs.(47), (^) we 
observe that the hamiltonian is of the form Hi = K + h, with for the specific 
case at hand a kinetic-energy term: 



^ = 2 s'^P^Pi 



1 



„2 , pI 

2 l^^+sin2( 



(55) 



and the potential (normalized for later convenience such that 2H = ujijOi): 



h{9,ip) = ho 



uj^p^ 



LOOJl 



ijj/uJl 



2pl 



1 — sin 0cos2((/9 — (/?*) 



(56) 



The corresponding lagrangean Li = K — h produces the Euler-Lagrange equa- 
tions 



Pe 



Pip 



„ n w^ sin 6* COS cos^ (cj9 — 09* ) 
sm 6* cos 6* c^'' + yr- ^ J 



d_ 

It 



(sin^ 6 ip) 



1 — sin^ 9cos'^{ip — 99*))^ ' 
w^ sin^ 6 sm{(p — ip^,) cos{ip — y?* 



(57) 



(1 — sin^ 9cos'^{ip — (p*)y 
These equations have solutions 

cos 6* = sin/3 sinwit, ^siv ~ ^*) = tg/3 coswit, 

with P a constant, implying the relation 

sm6 cos(99 — (p<t) = cos/?. 
Solving for the velocity (and taking into account the unit mass) 

Pe = V = — (Ji sin((^ — 9?^,), p^ = sin ^u''' = — cji cos/3 cos6', 



(58) 
(59) 
(60) 



in agreement with (pq). From these results we can compute Hamilton's principal 
function 



/■* 1 

Si{9,Lp,t)= / (irLi[6'(T),(/7(r)] = — arcctg^(tg6'sin(93-(y9*)) 
JO 2t 

This function indeed satisfies the Hamilton-Jacobi equations 



dSi 

do 



Pe. 



dSi 
dip 



Pip, 



with {pe-iPip) as given by eq.(^), and 



dSi 



OJOJl 



dSi 



+ 



1 fdSiy 



sm 



e\dipj 



hie, if). 



(61) 



(62) 



(63) 



Using the relation ctga = tg9sin{(p — ip^) = ctgwit, the equations ( |62D can 
be recast in the form 

Pi = LoViU = — ViT. (64) 

P* 

Hence T/p* can indeed be identified with Hamilton's principal function of this 
system. 

Repeating the arguments of sect. |2|, the action (21) for the Hamilton-Jacobi 
theory is now generalized to: 

A{p,Si-h) = Jdtjd^x^p (^^ + ^gid^,SiV,Si + h^ . (65) 

Reinterpretation of Si as a velocity potential for fluid flow: v = VS'i, leads 
back directly to the inhomogeneous Euler equation 

—l + vVvi = -Vih, (66) 



which for stationary flow becomes eq. (^q) . Variation of this action w.r.t. Si 
gives the equation of continuity for p, as before; note that in this action h plays 
the role of an external source for the density p. 

5 Quantum theory 

The quantum theory of a particle on a curved manifold is well-established. 
For the wave function to be well-defined and single- valued, the momenta must 
satisfy the Bohr-Sommerfeld quantization conditions 



® pidx^ 



2imh, (67) 



for any closed classical orbit T. For the free particle of unit mass on the unit 
sphere the left-hand side is 

pT 

v^dT = u?T = 2ttuj, (68) 

where T = 2ti /lo is the period of the orbit. Hence the quantization rule amounts 
to quantization of the rotation frequency (the angular momentum): uj = nfi. 
For the dual model, the same quantity takes the value 

- Vidx^ = V dr = — 2~;- = ujujiTi = 2Truj, (69) 

r Jo sin p 

and again uj = nfi. As the quantization conditions in the two dual models are 
the same, the duality can be preserved in the quantum theory. 

If this is to be true also in the fluid interpretation, the quantization con- 
ditions must be respected at that level as well. Now the first quantization 
condition for the integral ( |68| ) is interpreted in the fluid dynamical context as 
a quantization of the fluid momentum, cf. eq.(|3^). The second quantization 
condition (|69| ) has a twofold interpretation: first, according to eqs.(|4^), (p^) it 



quantizes the strength of the fluid sources and sinks in the model of free geodesic 
flow; the agreement between the two quantization conditions is then obvious: 
in order for the strength of the source/sink to satisfy a quantization condition, 
the amount of fluid transfered from one to the other must be quantized as well. 
In the context of the dual model however, the condition imposes the quan- 
tization of vorticity in the quantum fluid H. In the more general context of 
quantum models of fluids in geodesic flow on a compact two-dimensional surface 
and their duals described by the stream functions, this observation shows that 
duality at the quantum level requires the quantization of sources in one model 
to be directly related to the quantization of vorticity in the dual one. This 
situation closely parallels the relation between the quantization of monopole 
charge [^ and the quantization of the magnetic flux of fluxlines ||5| in three 
dimensions. 
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